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ABSTRACT 


The  combination  of  moving  target  indicator  (MTI)  and  coherent  in¬ 
tegration  filters  can  provide  velocity  filtering  with  advantages  that  neither 
technique  can  provide  alone.  These  techniques  are  evaluated  as  a  special 
case  of  a  transversal  filter  which  lends  itself  to  analysis  using  matrix  algebra. 
This  approach  gives  a  simple,  easily  understood  insight  into  the  effects  of 
weighting  the  input  radar  data.  The  resulting  matrix  equations  are  evaluated 
and  the  results  are  plotted  for  some  special  configurations.  The  computer 
programs  which  were  written  as  a  part  of  this  research  are  not  restricted 
to  these  configurations  nor  to  the  type  of  clutter  spectrum  that  is  assumed. 

It  is  shown  that  the  performance  of  coherent  integrators  in  a  coherent 
noise  (clutter)  environment  cannot  be  predicted  with  a  concise  analytical 
expression.  In  particular,  the  weighing  of  the  input  data  and  the  covariance 
matrix  of  the  interference  determine  the  integrator  performance. 

AUTHORIZATION 

NRL  Problem  R02-29 
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PERFORMANCE  OF  CASCADED  MTI  AND  COHERENT  INTEGRATION 
FILTERS  IN  A  CLUTTER  ENVIRONMENT 


I.  INTRODUCTION 

Recent  advances  in  velocity  filtering  for  radar  systems  have  been  made  possible  mainly 
through  the  use  of  digital  signal  processing.  The  detection  of  moving  targets  in  the  presence 
of  strong  returns  from  fixed  objects  (clutter)  has  classically  been  accomplished  with  moving 
target  indicators  (MTI)  or  with  pulse  doppler  radars  (Chap.  17  and  19,  Ref.  1).  The 
flexibility  of  digital  signal  processing  allows  improved  implementation  of  both  techniques, 
as  well  as  processing  configurations  which  make  use  of  the  advantages  of  both. 

MTI  is  normally  used  with  low-pulse-repetition-frequency  (low-PRF)  radars.  The  MTI 
filter  rejects  clutter  by  means  of  a  notch  in  its  passband  centered  on  the  clutter  doppler 
spectrum.  A  single  output  provides  moving  target  detection  over  the  remaining  doppler 
spectrum  between  the  PRF  harmonics.  Pulse  doppler  is  usually  associated  with  high-PRF 
radars.  A  contiguous  bank  of  narrow-band  filters  is  used  to  detect  moving  targets  outside 
of  the  clutter  spectrum. 

A  pulse  doppler  radar  can  usually  achieve  a  greater  improvement  in  the  signal-to-clutter 
ratio  than  can  an  MTI.  This  improvement  is  needed  to  contend  with  the  higher  clutter 
levels  caused  by  clutter  foldover  at  the  range  ambiguities  of  the  high-PRF  pulse  doppler 
radar.  The  increased  clutter  rejection  results  from  the  high  PRF  and  from  the  larger 
number  of  pulses  processed  by  the  pulse  doppler  radar. 

McAulay  (2)  formulated  the  MTI  problem  as  a  classical  detection  problem.  He 
determined  the  optimum  process  by  maximizing  the  resulting  likelihood  ratio  and  showed 
that  the  optimum  receiver  structure  could  be  interpreted  as  a  clutter  filter  in  cascade  with 
a  narrow-band  doppler  filter  bank.  This  arrangement  has  the  advantage  that  the  dynamic 
range  at  the  input  to  the  narrowband  processor  is  greatly  reduced  by  the  clutter  rejection 
of  the  clutter  filter. 

A  practical  approximation  to  this  processor  consists  of  a  conventional  MTI  filter 
cascaded  with  a  coherent  integrator  formed  by  a  contiguous  bank  of  narrow-band  filters. 

The  coherent  integrator  has  been  implemented  by  using  the  Fast  Fourier  Transform  (FFT) 
algorithm  (3,4).  The  processing  gain,  i.e.,  the  improvement  in  signal-to-clutter  ratio  that 
can  be  achieved  by  this  configuration,  is  considered  in  this  report. 

The  improvement  factor  that  can  be  achieved  using  an  MTI  alone,  and  then  a  coherent 
integrator  alone,  is  presented  first.  This  gives  a  basis  for  comparing  the  advantages  of 
cascading  the  two  filters.  A  Gaussian  .clutter  spectrum  is  assumed  in  computing  the 
improvement  factors.  Although  not  completely  accurate,  this  assumption  is  generally 
made  in  the  MTI  literature  and  it  has  been  found  to  give  a  reasonable  prediction  of  MTI 
performance. 
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In  general,  the  number  of  pulses  that  can  be  coherently  processed  is  limited  by  the 
number  of  pulses  that  are  transmitted  during  the  antenna  dwell  time.  This  number  can 
range  from  one  or  two  to  several  hundred,  depending  on  the  particular  radar  application. 
In  practice,  the  number  of  pulses  may  be  constrained  by  the  cost,  weight,  or  size  of  the 
processor  hardware.  In  this  report,  the  MTI  filters  are  considered  to  be  single,  double, 
or  triple  cancellers.  The  coherent  integrators  are  considered  to  be  8-  or  16-pulse  FFTs; 
the  weights  used  to  shape  the  FFT  filter  response  are  either  uniform  (all  weights  equal)  or 
Tschebyschev  (-25  dB  sidelobes).  These  choices  seem  to  be  typical  of  the  processors  that 
might  be  considered  for  many  low-PRF  radars.  However,  the  computer  programs  used  to 
evaluate  these  processors  are  neither  restricted  by  these  choices  nor  by  the  assumed 
Gaussian  clutter  spectrum. 


II.  GENERALIZED  PROCESSOR  MODEL 
Matrix  Equivalent 

Both  the  MTI  and  the  coherent  integrator  can  be  modeled  as  a  transversal  filter  of 
the  type  shown  in  Fig.  1.  A  tapped  delay  line  (or  multiple  delay  lines)  is  used  to  provide 
n+1  time  samples  of  the  input  ;c(t).  The  output  y(£)  consists  of  a  weighted  sum  of  these 
samples: 


y(f)  =  Wxxt  +  W2x(t  —  T)  +  ...  +  lFn+1  x(£  —  nT)  (1) 

where  T  is  the  time  between  samples.  x(t),  y(t),  and  the  weights  Wt  may  be  complex. 


■y(t) 


Fig.  1  —  Generalized  model  for  an  MTI  or  coherent  integrator  processor  consisting  of  a 

single  filter  and  a  single  output 


For  purposes  of  analysis,  it  is  convenient  to  write  Eq.  (1)  as  a  matrix  operation. 
Therefore,  define  W  and  X  as  column  vectors  with  elements  given  by 
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where  W,  is  the  ith  weight  shown  in  Fig.  1,  and 

JCj  =  *(t) 

x2  =  x(t-T) 

x3  =  x(t  —  2  T) 


xn+1  =  x(t-nT). 

With  these  matrices,  Eq.  (1)  can  be  written  as  the  product  of  the  transpose  of  W,  WT, 
times  X : 

y  =  WTX. 

The  output  power,  as  shown  in  Ref.  5,  is  given  by 

p  =  \y\2  =  yy*  =  yy* 

=  WtX(WtX)t  (2) 

=  WTXX^W 

where  use  has  been  made  of  the  fact  that 

yT  =  y. 

The  output  signal-to-clutter  ratio  can  be  computed  by  taking  the  ratio  of  Eq.  (2)  when 
x(t)  is  the  input  signal  to  Eq.  (2)  when  x(f)  is  the  input  clutter. 

Expected  Output  Clutter  Power 

The  expected  output  clutter  power  Pc  is  found  from  the  expected  value  of  Eq.  (2): 
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Pc  =  E  \WtXXtW*}  =WT  XXt  W*.  (3) 

Letting 

Mc  =  XX*T 

Then 


Pc  =  WT  Mc  W* 


(4) 


where  Mc,  the  covariance  matrix  of  the  clutter,  can  be  derived  from  the  Fourier  transform 
of  the  normalized  clutter  power  density  spectrum.  For  example,  using  a  Gaussian  clutter 
spectrum  yields 


Pdf)  = 


1 

e~2 


(5) 


The  mean  pc  implies  relative  motion  between  the  radar  platform  and  the  clutter.  The 
standard  deviation  oc  of  the  clutter  spectrum  is  a  measure  of  the  bandwidth  of  the  clutter 
spectrum. 


The  Fourier  transform  of  Eq.  (5)  is 

i P(t)  =  exp  (—27 r2oj?  r2  —  j  2ir  juc  r). 
Assuming  stationarity,  the  k,l  element  of  Mc  is 


Mc  (k,l)  =  exp[-27r2c;2  ( k  -  l)2Tz  -  j2tt  j uc  (k  -  Z)  T] .  (6) 


Expected  Output  Signal  Power 

Letting  the  input  x(t)  correspond  to  a  signal  input,  the  expected  output  signal  power 
Ps  is  found  in  the  same  way  as  above  for  the  clutter: 

Ps  =  WT  Ms  W*  (7) 

where  Ms  is  the  covariance  matrix  for  the  signal.  The  spectrum  of  a  target  is  affected  by 
such  things  as  the  target’s  characteristics,  transmitted  waveform,  radar  stability,  antenna 
scanning,  and  relative  velocity.  The  received  signal  waveform  r(f)  is  given  by 

r(t)  =  s{t)  eiWdt  (8) 

where  s(t)  includes  all  the  above  effects,  except  relative  velocity,  and  is  the  doppler 
frequency  corresponding  to  the  relative  velocity  between  radar  and  target.  Since  no  apriori 
knowledge  is  assumed  about  c od ,  it  is  given  a  uniform  probability  distribution  (0  <  cod  <  PRF) 

If  the  doppler  shift  is  the  dominant  effect,  i.e.,  the  bandwidth  of  the  other  effects 
are  small  relative  to  the  doppler  shifts  expected,  then  Eq.  (8)  is  approximately 
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•  C Jj  t 

r(t)  =eJ  d  .  (9) 

The  autocorrelation  of  Eq.  (9)  is 

i Hr)  =  e>“dT .  (10) 

The  processor  shown  in  Fig.  1  depicts  only  one  filter  and  one  output.  However,  in 
general,  multiple  independent  outputs  are  possible  by  selecting  appropriate  weight  vectors. 
In  particular,  the  FFT  weights  form  n+ 1  filters  with  n+1  outputs.  The  signal  gain  of  a 
particular  filter  then  will  be  the  expected  gain,  provided  that  ojd  is  such  that  a  signal 
exists  somewhere  between  the  crossover  points  of  adjacent  filters.  This  is  shown  in  Fig.  2 
for  an  8-pulse  FFT. 


Fig.  2  —  Signal  amplitude  vs  PRF  for  a  multiple-filter,  multiple-output 
integrator.  The  case  illustrated  here  is  for  an  8-pulse  FFT. 


From  Fig.  2,  the  gain  of  the  second  filter  is  derived  by  considering  a  target  doppler 
fd  that  can  occur  with  equal  probability  in  the  region 

f2  -  (B/ 2)  <  fd  <  f2  +  (B/2) 

where  f2  is  the  center  of  the  response  of  the  second  filter. 

The  gain  of  the  ith  filter  is  derived  in  the  same  way  by  considering  the  region 


fi  ~  (B/2)  <  fd  <  U  +  (B/2). 

The  probability  density  function  for  fd  associated  with  the  ith  filter  is 

'  1  IB,  for  (fi  -  B/2)  <fd<  (fi  +  B/2 

P(fd)  =  ' 

.0,  elsewhere. 

The  expected  value  of  Eq.  (10)  is 


E 


{m}  =  r 

J - c 


ei2v  f‘r  P{fd)  dfd  =  e12’  Si^(*gr). 


(11) 
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Assuming  stationarity,  the  k,l  element  of  the  signal  covariance  matrix  M$  is  given  by 


ms(k,l)  =  exp  [j 2tt  f{(k  -  l)  T] 


sin  [irBjk  -  l)  T ] 
irB(k  —  l)  T 


(12) 


Improvement  Factor 

The  processing  gain  of  any  arbitrary  transversal  filter  specified  by  a  weight  vector  IT 
can  be  computed  by  taking  the  ratio  of  the  signal  output  power  [Eq.  (7)]  to  the  clutter 
output  power  [Eq.  (4)] .  The  processing  gain,  or  improvement  factor  /,  is 


;  WT  Ms  W* 
WT  Mc  W*  ’ 


(13) 


The  matrix  Ms  is  generated  from  Eq.  (12),  and  Mc  is  generated  from  Eq.  (6). 


III.  MTI  FILTERS 

A  simplified  diagram  of  an  rc -stage  MTI  is  shown  in  Fig.  3.  This  diagram  is  equivalent 
to  the  model  in  Fig.  1  if  the  weights  are  given  by  the  binominal  coefficients  with  alternating 
signs,  i.e.. 


=  (-  1)*'  -  1  n_  ^  i  =  1,  2,  ....  n+1. 

Using  these  weights  in  Eq.  (7),  the  power  transfer  function  can  be  computed  by  setting 
the  bandwidth  B  equal  to  zero  in  Eq.  (12)  and  varying  from  zero  to  1/T,  where  1/T 
is  the  PRF.  This  corresponds  to  a  sine  wave  input  varying  from  zero  to  the  PRF.  The 
normalized  power  transfer  function  Ps,  as  defined  in  Eq.  (7),  is  found  for  n  =  1  to  n  =  7, 
and  the  results  are  as  shown  in  Fig.  4. 


Fig.  3  —  An  n -stage,  or  n -canceller,  MTI 
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Fig.  4  —  Normalized  power  transfer  function  Pg  for  an  MTI 
having  the  indicated  number  of  cancellers 


The  power  transfer  function  is  normalized  by  dividing  by  the  maximum  gain  WN  given  by 

n+ 1 

WN  -  ^  wf. 
i= 1 

The  maximum  gain  WN  is  listed  below  for  n  =  1  to  n  =  7: 


n 

wN 

1 

2 

2 

6 

3 

20 

4 

70 

5 

252 

6 

923 

7 

3432 

The  improvement  factor  for  an  n -stage  MTI  is  obtained  by  using  Eq.  (13).  The 
matrix  Ms  is  generated  by  letting  B  =  1/T  in  Eq.  (12),  i.e.,  the  target  doppler  may  occur 
at  any  frequency  from  zero  to  the  PRF  (PRF  =  1/T).  Also,  let  =  1/2 T,  i.e.,  the  center 
of  the  filter  is  at  one-half  the  PRF.  The  covariance  matrix  Mn  is  generated  using  Eq.  (6), 
with  the  average  clutter  doppler  nc  equal  to  zero.  The  results  are  plotted  in  Fig.  5,  where 
the  clutter  spectral  width  oc  ranges  from  10— 3  to  10— 1  times  the  PRF.  Curves  similar 
to  Fig.  5  have  been  derived  using  various  methods  and  published  in  several  places  (6,7). 
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CLUTTER  SPECTRAL  WIDTH/RRDRR  PRF 

Fig.  5  —  Improvement  factor  I  [given  by  Eq.  (13)]  for  an  MTI 
having  the  indicated  number  of  cancellers  (filters) 


IV.  COHERENT  INTEGRATOR  (FFT) 
Uniform  Weights 


The  discrete  Fourier  transform  can  be  characterized  by 


n  =  l 


(14) 


where  xn  is  the  input  time  sample,  and  the  yk  are  the  output  spectral  components. 
Therefore,  Eq.  (14)  describes  N  filters  with  outputs  y1(  y2>  •••»  as  sketched  in  Fig.  2. 
The  effective  weights  of  the  ftth  filter  are 


— y  2  tt  ~  D 

Wkn=e  N  ,  n  =  1,  2,  ....  N.  (15) 

These  filters  have  a  (sin  x)/x  shape  and  are  translated  in  frequency  by  the  index  ft. 
When  N  =  16,  the  filter  responses  for  ft  =  1  and  ft  =  2  are  as  shown  in  Figs.  6(a)  and  6(b), 
respectively.  The  normalization  used  for  these  curves  is 

N 

WN=^\Wkn\*  =  256 

n  =  i 


for  N  =  16. 


FILTER  GRIN/MRXo  FILTER  GRIN  (OBI  FILTER  GRIN/MRX„  FILTER  GRIN  (DB) 
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(a) 


(b) 

Fig.  6  —  Normalized  power  transfer  function  (filter  response)  for  (a) 
filter  no.  1  and  (b)  filter  no.  2  of  a  16-pulse  (16-filter)  MTI  integrator 
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Improvement  Factor  —  The  improvement  factor  for  the  coherent  integration  is 
computed  the  same  way  as  for  the  MTI  discussed  in  Section  III.  Equation  (13)  is  used 
with  the  weight  vector  W  generated  by  Eq.  (15).  The  signal  covariance  matrix  Ms  is 
generated  using  Eq.  (12),  with  the  bandwidth  B  given  by 


where  N  is  the  number  of  pulses  integrated.  The  center  of  the  ith  filter,  /),  is  given  by 

fi  =  -4r,  i  =  1,  2,  ...,  N. 

The  improvement  factor  is  shown  in  Fig.  7(a)  for  an  8-pulse  integrator  and  in  Fig.  7(b) 
for  a  16-pulse  integrator.  The  gains  of  the  filters  are  averaged  and  shown  by  a  dotted 
curve  in  each  figure.  The  average  gain  corresponds  to  the  expected  improvement  factor 
for  a  target  with  equal  probability  of  occurring  in  any  filter.  This  average  improvement 
factor  should  be  compared  with  the  MTI  improvement  factor  shown  in  Fig.  5,  which  was 
computed  for  a  target  with  equal  probability  of  occurring  at  any  doppler. 

The  irregularities  in  Fig.  7  are  caused  by  the  interaction  between  the  clutter  spectrum 
and  the  filter  transfer  function.  Note  that  the  width  of  the  sidelobes  of  the  filters,  as 
shown  in  Fig.  6,  is  given  by  PRF /N,  i.e.,  by  0.125  X  PRF  for  an  8-pulse  integrator  and  by 
0.0625  X  PRF  for  a  16-pulse  integrator.  At  the  left  side  of  Fig.  7  the  clutter  spectrum  is 
narrow  as  compared  to  the  sidelobe  width  so  that  more  of  the  clutter  energy  coincides 
with  the  null  between  sidelobes.  The  improvement  factor  is  reduced  as  the  clutter  width 
is  increased  and  more  energy  coincides  with  a  peak  in  the  sidelobes.  This  effect  continues 
until  the  clutter  width  becomes  wide  enough  for  an  appreciable  amount  of  clutter  energy 
to  coincide  with  the  next  sidelobe  null.  In  this  region,  less  clutter  energy  is  passed,  and 
the  improvement  factor  decreases  at  a  slower  rate  as  the  clutter  spectral  width  is  increased. 
This  effect  continues  until  an  appreciable  amount  of  clutter  energy  coincides  with  the 
next  peak  in  the  sidelobes. 

For  narrow  clutter  spectral  widths,  the  improvement  factors  for  the  8-  and  16-pulse 
integrators  are  approximately  the  same,  as  can  be  seen  in  Figs.  7  (a)  and  (b).  Intuition 
would  seem  to  indicate  that  the  integration  of  more  pulses  should  result  in  more  gain  in 
the  signal-to-clutter  ratio.  When  the  interference  consists  of  “white”  noise,  the  improve¬ 
ment  I  in  signal-to-noise  ratio  is  given  by 

/  =  10  log  N 

which,  as  expected,  does  give  a  larger  improvement  as  N  is  increased. 

For  correlated  (or  “colored”)  noise,  the  gain  that  can  be  realized  depends  on  the 
shape  of  the  noise  spectrum  and  the  shape  of  the  filter  transfer  function.  Except  for 
filter  number  one,  all  the  filters  generated  by  both  the  8-  and  16-pulse  integrators  have 
sidelobe  nulls  at  zero  doppler.  The  width  of  the  null  for  an  8-pulse  integrator  is  twice 
as  wide  as  for  a  16-pulse  integrator.  Apparently,  the  additional  clutter  rejection  caused 
by  the  wider  null  of  the  8-pulse  integrator  offsets  the  additional  integration  gain  of  the 
16-pulse  integrator. 


CLUTTER  SPECTRAL  WIBTH/RRDAR  PRF 

(b) 

Fig.  7  —  Improvement  factor  I  [given  by  Eq.  (13)]  for  the  indicated 
filter  number  N  in  (a)  an  8-pulse  and  (b)  a  16-pulse  coherent  integrator 
with  uniform  weighting.  The  average  improvement  for  all  filters  is 
indicated  by  the  dotted  curve. 
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For  wide  clutter  spectral  widths,  this  effect  disappears  and  the  16-pulse  integrator 
achieves  more  gain.  The  wider  spectral  width  is  more  characteristic  of  “white”  noise. 

Notice  that  filter  number  one  has  its  mainlobe  at  zero  doppler,  as  shown  in  Fig.  6(a). 
The  clutter  spectrum  is  also  centered  at  zero  doppler.  For  these  reasons,  the  improvement 
factor  for  filter  number  one,  as  shown  in  Fig.  7,  increases  as  the  clutter  spectral  width 
increases  since  this  results  in  more  clutter  energy  outside  the  filter  mainlobe  region. 


Tschebyschev  Weights 

The  sidelobes  levels  of  the  filters  shown  in  Fig.  6  can  be  reduced  by  using  a  data 
“window”  function  (Ref.  8).  This  corresponds  to  weighting  the  input  time  samples  xn 
of  Eq.  (14).  Let  xn  -*  an  xn  in  Eq.  (14)  where  an  is  the  weight  vector  for  a  desired 
filter  shape.  Then,  from  Eq.  (15),  the  total  effective  weights  of  the  feth  filter  become 

—j2n  il^iy.j-1) 

Wkn=ane  N  ,  n  =  1,  2,  ....  N.  (16) 

When  the  weights  an  correspond  to  the  familiar  Tschebyschev  weights  for  a  —25  dB 
sidelobe  level,  a  16-pulse  coherent  integrator  results  in  16  filters,  two  of  which  are  shown 
in  Figs.  8(a)  and  (b). 

Improvement  Factor  —  The  improvement  factors  for  8-  and  16-pulse  integrators 
with  —25  dB  sidelobes  are  shown  in  Figs.  9(a)  and  9(b),  respectively.  Notice  that  at 
narrow  clutter  spectral  widths,  the  average  improvement  factor  of  the  8-pulse  integrator 
is  greater  than  the  average  improvement  of  the  16-pulse  integrator.  Also,  notice  the 
irregular  ordering  of  the  filters  in  Fig.  9(b)  for  narrow  clutter  spectral  widths.  Those 
filters  whose  center  frequencies  are  farther  from  the  clutter  spectrum  center  (zero  doppler) 
might  be  expected  to  have  a  higher  improvement  factor  than  filters  nearer  to  the  clutter. 
However,  for  a  narrow  clutter  spectrum,  this  does  not  happen.  These  peculiarities  are  a 
result  of  the  irregular  sidelobe  and  null  widths  produced  by  the  Tschebyschev  weighting. 

A  filter  that  is  closer  to  the  clutter  may  have  a  wider  null  near  zero  doppler  and,  therefore, 
realize  a  higher  improvement  factor. 

An  examination  of  Figs.  7  and  9  reveals  a  certain  degree  of  unpredictability  for 
coherent  integration  filters  in  a  clutter  (or  “colored”  noise)  environment.  This  unpredict¬ 
ability  could  be  reduced  by  preceding  the  integration  filters  with  a  “clutter”  filter,  such 
as  an  MTI  canceller,  to  reduce  the  coherency  of  the  clutter  input  to  the  integration  filters. 
This  configuration  will  be  examined  in  the  following  section. 


V.  CASCADED  MTI  AND  COHERENT  INTEGRATOR 
Optimum  Cascaded  Filters 

The  maximum  improvement  in  signal-to-clutter  ratio  could  be  achieved  by  following 
a  “prewhitening”  filter  with  a  filter  matched  to  the  target  return.  Since  the  target  doppler 
spectrum  considered  by  this  report  is  a  single  line  component,  the  matched  filter  cannot 


DOPPLER  FREQUENCE /RADAR  PRF 


(b) 

Fig.  8  —  Normalized  power  transfer  function  (filter  response)  for  (a) 
filter  no.  1  and  (b)  filter  no.  2  in  a  16-pulse  coherent  integrator.  The 
input  time  samples  xn  are  Tschebyschev  weighted  for  a  —25  dB  side- 
lobe  level.  (Compare  with  Fig.  6.) 


IMPROVEMENT  FRCTOR  (DB) 
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(a) 


(b) 

Fig.  9  —  Improvement  factor  I  [given  by  Eq.  (13)]  for  (a)  an  8-pulse 
and  (b)  a  16-pulse  coherent  integrator  with  Tschebyschev  weighting. 
The  average  improvement  for  all  filters  is  indicated  by  the  dotted 
curve. 
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be  realized  by  a  processor  with  a  finite  processing  time.  However,  the  coherent  integration 
filters  considered  in  the  previous  section  are  an  approximation  to  a  matched  filter.  The 
approximation  is  limited  only  by  the  processing  time  (or  the  number  of  pulses  integrated). 

To  completely  “prewhiten”  the  clutter  prior  to  coherent  integration,  the  prewhitening 
filter  must  have  a  processing  time  equal  to  the  processing  time  of  the  coherent  integrator, 
i.e.  each  filter  must  process  the  same  number  of  pulses.  Although  this  approach  gives  the 
classical  matched  filter  in  “colored”  noise,  and  therefore  gives  the  maximum  improvement 
in  signal-to-clutter  ratio,  it  is  worthwhile  to  consider  substituting  for  the  ideal  prewhitening 
filter  a  more  simple  filter  with  a  shorter  processing  time.  Any  clutter  notch  filter  that 
precedes  the  integration  filters  reduces  the  dynamic  range  of  the  input  to  the  integration 
filters  by  an  amount  equal  to  the  clutter  rejection  realized  by  the  clutter  notch.  This 
requires  less  storage  for  the  integration  filters. 

In  this  report,  the  three  clutter  notch  filters  considered  are  MTIs  consisting  of  (a)  a 
single  canceller,  (b)  a  double  canceller,  and  (c)  a  triple  canceller.  The  integration  filters 
considered  are  an 8- and  a  16-pulse  FFT.  Therefore,  the  processing  time  of  the  clutter 
filter  is  less  than  the  processing  time  of  the  integration  filter. 

Some  care  is  required  in  the  analysis  of  a  cascaded  combination  of  two  sampled  data 
filters  with  different  processing  times.  Theoretically,  the  order  of  the  filters  can  be  inter¬ 
changed  with  the  same  results  since  both  filters  are  linear.  For  analysis  it  is  convenient 
to  consider  the  filter  with  the  shorter  processing  time  first,  and  then  the  filter  with  the 
longer  processing  time.  Also,  for  practical  reasons,  the  MTI  should  precede  the  integration 
filter  to  reduce  the  dynamic  range  that  the  integration  filter  has  to  handle. 


Matrix  Equations  For  Cascaded  Filters 

To  analyze  the  cascaded  filters  in  terms  of  the  matrix  operation  of  the  preceding 
sections,  consider  the  filter  arrangement  shown  in  Fig.  10.  The  first  filter  is  characterized 
by  the  number  M  of  pulses  processed  and  by  the  weight  vector  A .  The  second  filter  is 
characterized  by  the  number  N  of  pulses  it  processes  and  by  its  weight  vector  B.  The 
number  M  is  always  less  then  N. 

The  output  z(t)  of  these  filters  can  be  expressed  in  terms  of  the  input  y(t)  to  the 
second  filter  by  means  of  the  matrix  equation 

z=BtY  (17) 


Vl~ 

*r 

y  2 

b2 

Y  = 

. 

and  B  = 

• 

yN_ 

_bN_ 

where 
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z(f) 


Fig.  10  —  Cascaded  filters  with  different  processing  times.  M  and 
N  are  the  number  of  pulses  processed  ( M  <  N),  and  A  and  B  are 
weight  vectors. 


with  Bt  being  the  transpose  of  the  weight  vector  B.  The  components  of  Y  can  be  related 
to  the  input  x(t)  by 


Vl  "  ai*i  +  a2*2  +  -  +  aM  XM 
y2  =  a\x2  +  a2x3  +  —  +  aM  XM+ 1 


yN  -  atxN  +  a2xN+1  +  ...  +  aM  *m+N-1- 


Let 


*1 

x2 


X  = 


and  A1  = 


xM+N—1 


aM 


0 

0 


N— 1  zeros 


Then 


Vl  =(A!)tX. 


In  the  same  way 


O' 2  “  (^2  )r^» 
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where 


An  - 


aM 

0 

0 


1  zero 


N— 2  zeros 


and 

y n  =  ( an)t  x 

where 


0 

0 


N— 1  zeros 


Combining  all  of  these  matrix  relations  gives 

Y  =  AtX 


(19) 


where  X  is  given  by  Eq.  (18)  and 

A  =  [AiA2—An'\  . 

Substituting  Eq.  (19)  into  Eq.  (17),  the  matrix  relationship  between  input  and  output 
becomes 
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since  the  total  effective  weight  vector  is  given  by 

W  =  AB.  (21) 

Therefore  this  cascaded  combination  of  filters  can  be  evaluated  in  the  same  way  that 
was  outlined  in  section  II  and  used  in  section  III  for  the  MTI  filter  and  in  section  IV  for 
the  integration  filter.  The  only  difference  is  that  the  weight  vector  IV  is  defined  by  Eq.  (21) 
and  generated  as  shown  above. 


Transfer  Function 

The  power  transfer  function  of  the  cascaded  filters  can  be  computed  by  means  of 
Eq.  (7)  and  the  weights  given  by  Eq.  (21).  The  signal  covariance  matrix  Ms  is  generated 
from  Eq.  (12)  by  setting  the  bandwidth  B  equal  to  zero  and  varying  the  center  frequency 
fi  from  zero  to  1  IT.  The  weight  vector  IV  is  computed  from  Eq.  (21),  with  the  matrix  A 
generated  from  the  weight  vector  A  given  by 


al 

a2 


If  the  first  filter  is  an  MTI  consisting  of  M  —  1  binomial  cancellers,  then 

ai  =  (~  1 Y  ~  1  (Yri),  «'  =  1»  2,  ...,  M  (22) 

as  shown  in  section  III.  If  the  second  filter  is  an  FFT  with  uniform  weights,  then  the 
weight  vector  B  can  be  generated  from  Eq.  (15). 


Using  this  procedure,  the  power  transfer  functions  generated  by  a  triple-canceller 
MTI  cascaded  with  a  16-pulse  coherent  integrator  are  shown  in  Figs.  ll(a)-(i)  for  nine 
of  the  sixteen  filters.  Filter  no.  10  is  the  mirrow  image  of  filter  no.  8;  filter  no.  11  is 
the  mirrow  image  of  filter  no.  7,  etc.  Uniform  integrator  weights  are  used.  When  25-dB 
Tschebyschev  weights  are  used,  the  resulting  transfer  functions  are  shown  in  Figs.  12(a)-(i). 
The  filter  patterns  are  analogous  to  the  antenna  patterns  obtained  with  a  16-element  linear 
array  with  elements  equally  spaced  at  one-half  wavelength.  The  element  pattern  is 
analogous  to  the  MTI  filter  transfer  function.  All  filter  transfer  functions  are  normalized 
to  the  maximum  gain  of  the  center  filter(filter  no.  (N/ 2)  +  1,  where  N  is  the  number  of 
pulses  integrated). 


Improvement  Factor 

The  improvement  in  signal-to-clutter  ratio  is  shown  in  Figs.  13(a)  and  (b)  for  a  single 
canceller  cascaded  with  an  8-pulse  integrator.  With  25-dB  Tschebyschev  weights  [Fig.  13(b)] 
there  is  a  larger  variation  in  gain  between  filters.  However,  the  average  gain  is  greater  by  5  to 
6  dB  than  for  uniform  weights  [Fig.  13(a)] .  The  integrator  with  uniform  weights  has  more 
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gain  in  the  filters  near  the  clutter  spectrum,  i.e.,  filters  no.  2  and  8.  The  MTI  null  at  zero 
doppler  overrides  the  gain  that  was  achieved  by  filter  no.  1  without  MTI. 

Figures  14(a)  and  (b)  show  the  improvement  factors  for  a  single  canceller  cascaded  with 
a  16-pulse  integrator.  Again,  with  25-dB  Tschebyschev  weights  [Fig.  14(b)]  there  is  a 
larger  variation  in  gain  between  filters.  The  filters  near  the  clutter  spectrum  again  have 
less  gain  than  the  same  filters  with  uniform  weights  [Fig.  14(a)] .  However,  the  average 
filter  gain  is  slightly  higher  (about  3  dB),  for  narrow  clutter  spectral  widths  (0.001  X  PRF) 
using  uniform  weights,  but  slightly  higher  (about  2  dB)  for  wide  clutter  widths  (0.1  X  PRF) 
using  25-dB  weights.  Comparing  Figs.  13  and  14  shows  that  the  average  improvement 
factor  is  still  less  for  a  16-pulse  integrator  than  for  an  8-pulse  integrator  for  narrow  clutter. 
This  effect  was  noticed  for  the  8-  and  16-pulse  integrators  without  MTI. 

Figures  15(a)  and  (b)  show  the  improvement  factors  for  a  double  canceller  cascaded 
with  an  8-pulse  integrator.  The  characteristics  of  these  curves  are  similar  to  Fig.  13  for 
a  single  canceller  and  an  8-pulse  integrator.  With  25-dB  weights  [Fig.  15(b)]  there  is  a 
larger  variation  in  filter  gains,  but  also  a  larger  (again  5  to  6  dB)  average  filter  gain.  The 
additional  improvement  factor  realized  by  preceding  an  8-pulse  integrator  with  a  double 
canceller  (Fig.  15)  instead  of  with  a  single  canceller  (Fig.  13)  varies  from  42  dB  for  narrow 
clutter  widths  (0.001  X  PRF)  to  4  dB  for  wide  clutter  widths  (0.1  X  PRF)  using  either 
uniform  or  25-dB  weights. 

Figures  16(a)  and  (b)  show  the  improvement  factor  for  a  double  canceller  cascaded 
with  a  16-pulse  integrator.  Comparing  uniform  weights  [Fig.  16(a)]  with  25-dB  weights 
[Fig.  16(b)] ,  the  difference  in  the  variation  in  filter  gain  is  reduced,  although  the  filters 
near  the  clutter  spectrum  still  have  less  gain  with  25-dB  weights  for  narrow  clutter.  For 
narrow  clutter  widths,  uniform  weights  give  about  2  dB  more  average  gain.  For  wide 
clutter  widths,  25-dB  weights  give  about  2  dB  more  average  gain.  The  additional  improve¬ 
ment  factor  realized  by  preceding  a  16-pulse  integrator  with  a  double  canceller  (Fig.  16) 
instead  of  with  a  single  canceller  (Fig.  14)  varies  from  about  42  dB  for  narrow  clutter  to 
about  4  dB  for  wide  clutter  using  either  uniform  or  25-dB  weights.  This  is  the  same 
improvement  that  was  achieved  by  preceding  the  8-pulse  integrator  with  a  double  instead 
of  with  a  single  canceller. 

Figures  17(a)  and  (b)  show  the  improvement  factor  for  a  triple  canceller  cascaded 
with  an  8-pulse  integrator.  Comparing  uniform  weights  [Fig.  17(a)]  with  25-dB  weights 
[Fig.  17(b)] ,  the  25-dB  weights  still  cause  a  wider  variation  in  filter  gains,  with  the  filters 
near  the  clutter  having  appreciably  less  gain  than  the  same  filters  with  uniform  weights. 

The  average  gain  is  always  greater  with  25-dB  weights,  varying  from  about  6  dB  with 
narrow  clutter  to  about  4  dB  for  wide  clutter.  The  additional  improvement  factor  realized 
by  preceding  an  8-pulse  integrator  with  a  triple  canceller  (Fig.  17)  instead  of  with  a  double 
canceller  (Fig.  15)  varies  from  more  than  30  dB  for  narrow  clutter  to  about  3  dB  for  wide 
clutter  using  uniform  weights,  and  varies  from  more  than  30  dB  to  about  2  dB  using  25-dB 
weights. 

Figures  18(a)  and  (b)  show  the  improvement  factor  for  a  triple  canceller  cascaded 
with  a  16-pulse  integrator.  Again,  the  25-dB  weights  cause  a  wider  variation  in  filter  gains. 
The  average  filter  gain  is  always  about  2  dB  greater  for  25-dB  weights.  The  additional 
improvement  realized  by  preceding  the  16-pulse  integrator  with  a  triple  canceller  (Fig.  18) 
instead  of  with  a  double  canceller  (Fig.  16)  varies  from  more  than  30  dB  for  narrow 
clutter  to  about  3  dB  for  wide  clutter.  This  is  the  same  increase  that  was  achieved  by 
preceding  the  8-pulse  integrator  with  a  triple  canceller  instead  of  with  a  double  canceller. 
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OOPPLER  FREQUENCY /RflDflR  PRF  OOPPLER  FREQUENCY/RAOAR  PRF 


(c)  (d) 

Fig.  11  —  Normalized  power  transfer  function  for  a  triple-canceller  MTI  cascaded  with  a  16-pulse 
coherent  integrator.  The  transfer  functions  for  nine  of  the  sixteen  filters  (filters  no.  1  through 
9)  are  shown  in  the  graphs  labeled  (a)  through  (i).  The  integrator  filters  have  been  uniformly 
weighted. 
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(C)  (d) 

Fig.  12  —  Normalized  power  transfer  function  for  a  triple -canceller  MTI  cascaded  with  a  16-pulse 
coherent  integrator.  The  transfer  functions  for  nine  of  the  sixteen  filters  (filters  no.  1  through 
9)  are  shown  in  the  graphs  labeled  (a)  through  (i).  The  integrator  filters  have  been  Tschebyschev 
weighted  for  a  -25  dB  sidelobe  level  (Compare  with  Fig.  11.) 


'IlTER  GAIN/MAX.  KILTER  GAINIGBI  filter  grin/max.  filter  ghiniqbi 


NRL  REPORT  7533 


23 


IMPROVEMENT  FACTORto  IMPROVEMENT  FACTOR  (DB) 


G.  A.  ANDREWS,  JR. 


(a) 


(b) 


Fig.  13  —  Improvement  factor  I  [given  by  Eq.  (13)]  for  a  single -canceller 
MTI  cascaded  with  an  8-pulse  integrator.  The  integrator  filters  have  been 
(a)  uniformly  weighted  and  (b)  25-dB  Tschebyschev  weighted.  The 
average  improvement  for  all  filters  is  indicated  by  the  dotted  curve. 
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(a) 


CLUTTER  SPECTRAL  HIDTH/RADAR  PRF 
(b) 

Fig.  14  —  Improvement  factor  I  [given  by  Eq.  (13)]  fora  single-canceller 
MTI  cascaded  with  a  16-pulse  integrator.  The  integrator  filters  have  been 
(a)  uniformly  weighted  and  (b)  25-dB  Tschebyschev  weighted.  The  aver¬ 
age  improvement  for  all  filters  is  indicated  by  the  dotted  curve. 
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(a) 


CLUTTER  SPECTRAL  WIDTH/RADAR  PRF 


(b) 

Fig.  15  —  Improvement  factor  I  [given  by  Eq.  (13)]  for  a  double -canceller 
MTI  cascaded  with  an  8-pulse  integrator.  The  integrator  filters  have  been 
(a)  uniformly  weighted  and  (b)  25-dB  Tschebyschev  weighted.  The  aver¬ 
age  improvement  for  all  filters  is  indicated  by  the  dotted  curve.  (Com¬ 
pare  with  Fig.  13.) 
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CLUTTER  SPECTRAL  WIDTH/RRORR  PRF 
(b) 


Fig.  16  —  Improvement  factor  I  [given  by  Eq.  (13)]  fora  double -canceller 
MTI  cascaded  with  a  16-pulse  integrator.  The  integrator  filters  have  been 
(a)  uniformly  weighted  and  (b)  25-dB  Tschebyschev  weighted.  The  aver¬ 
age  improvement  for  all  filters  is  indicated  by  the  dotted  curve.  (Com¬ 
pare  with  Fig.  14.) 
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(a) 


(b) 

Fig.  17  —  Improvement  factor  I  [given  by  Eq.  (13)]  for  a  triple-canceller 
MTI  cascaded  with  an  8-pulse  integrator.  The  integrator  filters  have  been 
(a)  uniformly  weighted  and  (b)  25-dB  Tschebyschev  weighted.  The  aver¬ 
age  improvement  for  all  filters  is  indicated  by  the  dotted  curve.  (Com- 


(b) 

Fig.  18  —  Improvement  factor  I  [given  by  Eq.  (13)]  for  a  triple-canceller 
MTI  cascaded  with  a  16-pulse  integrator.  The  integrator  filters  have  been 
(a)  uniformly  weighted  and  (b)  25-dB  Tschebyschev  weighted.  The  aver¬ 
age  improvement  for  all  filters  is  indicated  by  the  dotted  curve.  (Com¬ 
pare  with  Fig.  16.) 
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VI.  CONCLUSIONS 

Coherent  integration  is  a  technique  that  has  been  used  in  radar  systems  to  improve 
detection  in  a  “white”  noise  (i.e.,  receiver  noise)  environment  and  also  to  reject  clutter 
in  high-PRF  radars.  More  recently  this  technique  has  been  considered  for  clutter  rejection 
in  low-PRF  radars  where  the  ratio  of  the  clutter  spectrum  width  to  the  PRF  is  much 
larger  than  for  high-PRF  systems.  Less  analysis  is  available  for  this  application  to  aid  in 
understanding  and  evaluating  the  performance  of  coherent  integration. 

This  report  has  considered  several  particular  doppler  filters  which  consist  of  cascaded 
MTI  and  coherent  integration  filters.  These  filters  are  considered  typical  of  the  filters  that 
would  be  applicable  to  a  low-PRF  radar  which  has  a  relatively  small  number  of  hits  per 
dwell  time.  The  computed  improvement  factors  presented  in  Figs.  7,  9,  and  13-18  show 
that  coherent  integration  filters  perform  in  an  irregular,  almost  unpredictable,  way  then 
their  input  is  “colored”  noise  (or  clutter).  The  performance  of  coherent  integrators  against 
“white”  noise  is  well  established  and  easily  predicted.  Against  “clutter”,  the  performance 
is  highly  dependent  upon  the  shape  of  the  clutter  spectrum  and  the. shape  (or  weights)  of 
the  coherent  integrator  transfer  functions.  In  particular,  it  is  not  necessarily  true  that 
more  improvement  is  achieved  by  integrating  more  pulses  for  a  given  weight  vector.  On 
the  other  hand,  it  is  theoretically  true  that  weights  exist  which  will  allow  more  gain  to 
be  achieved  when  more  pulses  are  processed.  The  computation  of  weight  vectors  that 
will  give  the  maximum  improvement  in  signal-to -clutter  ratio  under  various  optimization 
criteria  will  be  presented  in  a  later  report. 

The  optimum  weights  obviously  depend  on  the  optimization  criteria  as  well  as  on 
the  shape  of  the  interference  spectrum.  Of  the  integration  filter  weights  discussed  in  the 
report,  either  uniform  or  25-dB  Tschebyschev  weights  could  be  selected  as  the  “better” 
weights  depending  on  the  criterion  used.  This  is  pointed  out  by  the  following  specific 
results: 

a.  25-dB  Tschebyschev  weights  give  more  average  gain  than  uniform  weights. 

b.  Uniform  weights  generally  give  more  gain  in  the  filters  near  the  clutter  spectrum. 

c.  Uniform  weights  result  in  less  variation  in  gain  from  filter  to  filter. 

d.  For  very  narrow  clutter  spectral  widths,  the  integration  of  more  pulses  does  not 
necessarily  lead  to  more  average  gain. 

These  conclusions  apply  regardless  of  whether  or  not  an  MTI  precedes  the  integrator. 
When  an  MTI  is  cascaded  with  the  integrator,  the  following  general  conclusions  can  be 
made: 


a.  The  additional  average  gain  achieved  by  cascading  an  MTI  with  any  particular 
integrator  is  only  very  slightly  affected  by  the  type  of  integrator.  That  is,  cascading  a 
single  canceller  with  an  8-pulse  integrator  results  in  about  the  same  additional  gain  above 
an  8-pulse  integrator  alone  as  the  additional  gain  above  a  16-pulse  integrator  achieved  by 
cascading  a  single  canceller  with  a  16-pulse  integrator.  [Compare  Fig.  7(a)  with  Fig.  13(a) 
and  Fig.  7(b)  with  14(a)] .  An  equivalent  comparison  can  be  made  for  double  and  triple 
cancellers. 
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b.  The  average  gain  achieved  by  cascading  an  integrator  with  any  particular  MTI  is 
affected  to  some  degree  by  the  type  of  MTI.  This  can  be  seen  by  comparing  the  improve¬ 
ment  factors  for  the  MTI,  shown  in  Fig.  5,  with  those  in  Figs.  13-18  for  the  MTI  cascaded 
with  integrators. 

It  might  be  concluded  that  the  chief  advantage  of  the  coherent  integrator  lies  not 
with  its  average  performance  in  a  clutter-dominated  environment  but  with: 

a.  Increased  gain  for  targets  near  the  clutter  spectrum.  It  narrows  the  “blind” 
velocity  regions  and  could  increase  the  detection  of  low-speed  targets. 

b.  Increased  gain  in  a  “white”  (i.e.,  receiver)  noise  environment.  This  also  applies 
to  other  forms  of  wideband  interference  such  as  clutter  received  through  the  antenna 
sidelobes  when  the  radar  is  on  a  moving  platform.  MTI  alone  has  no  capability  against 
this  type  of  interference. 

c.  Improved  velocity  measurement  of  targets  for  tracking  or  identification. 

These  factors  must  be  considered  in  any  tradeoff  between  cancellation  and  integra¬ 
tion  for  a  particular  radar  application  and  for  a  particular  interference  environment. 
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